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SIGNIFICANCE  AND  EXPLANATION 


Polynomial  interpolation  to  a univariate  function  is  a thoroughly  studied 
subject  and  when  properly  executed  provides  a numerically  efficient  means  of 
approximation . 

Very  little  seems  to  be  known  about  polynomial  interpolation  of  multi- 
variate functions.  However,  Kergin  recently  established  the  existence  and 
uniqueness  of  a natural  extension  of  univariate  interpolation  to  a multivariate 
setting.  In  this  paper  we  provide  a formula  for  Kergin  interpolation.  This 
formula  is  based  on  the  Newton  form  for  univariate  polynomial  interpolation. 

The  error  in  approximating  by  Kergin  interpolation  is  also  obtained  in  a 
convenient  form  which  allows  us  to  assess  the  quality  of  this  scheme.  In 
particular , we  establish  that  Kergin  interpolation  converges  for  analytic 
functions  of  several  variables.  

There  seems  to  be  a close  connection  between  Kergin  interpolation  and  the 
multidimensional  B- spline  introduced  by  de  Boor.  In  univariate  approximation 
B-splines  play  a central  role  in  both  the  theory  and  numerical  application  of 
spline  functions.  An  extremely  useful  numerically  stable  recurrence  formula 
for  univariate  B-splines  is  known  and  widely  used.  He  provide  here  a similar 
formula  for  multivariate  B-splines. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive  sui 
lies  with  MRC,  and  not  with  the  author  of  this  report. 
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A CONSTRUCTIVE  APPROACH  TO  KEPGIN  INTERPOLATION  IN  R 
Charles  A.  Micchelli 

In  his  recent  doctoral  thesis  P.  Kergin  [3]  introduced  a novel  method  of  interpolating 
smooth  functions  of  several  variables.  In  this  note,  we  give  an  explicit  representation 
for  Kargin's  interpolation  scheme. 

He  let  w (Rk)  denote  the  space  of  polynomials  of  total  degree  < n and  C^(Rk) 
n — 

be  all  functions  with  n continuous  derivatives  on  r\ 

Kergin  proves  the  following 

Theorem.  Given  any  x°,...,xn  e Rk  , not  necessarily  distinct,  there  is  a unique  linear 
map  P : C^  -*■  with  the  property  that  for  every  f « C*n',  every  homogeneous  differ- 
ential operator  q of  order  t , 0 <_  l <_  n,  and  every  subset  J c_  {0,1,...  ,n>  with  card 
J ■ Z+l  there  exists  a point  x in  the  convex  hull  {x^|j  e J}  such  that  qPf  (x)  - qf(x). 

Pf  necessarily  interpolates  f at  x , ...,xn  as  may  be  seen  by  choosing  i-  0. 

If  a point  is  repeated  with  some  multiplicity  then  all  the  derivatives  of  Pf  agree  with 


f at  that  point  with  the  same  multiplicity.  In  particular,  when  all  the  points  coincide 

x^  ■ x°  then  Pf  is  the  Taylor  polynomial  of  f at  x°  , 

Pf(x)  - f(x°)  + D . f(x°)  + ...  + \ Dn  _f<x°). 

0 nl  0 

x-x  x-x 

In  this  formula,  0 Qf  is  the  directional  derivative  of  f in  the  direction  x-x°. 
x-x 

For  the  general  form  of  Kergin' s map  we  use  the  notation 


/ f(v„x°  + v.x1  + ...  tv  xn)dv,  ...dv 
0 1 n 1 n 

sn 


Sn  - {v  - (vQ,...,vn)  |vj  ^0,  j v - 1}  . 

Me  will  show  that  Kergin  *s  map  is  given  by 

(1)  Pf(x)  - f(x°)  + / D f + ...  + ' / 


, 0 1,  x-x 

(x  ,x  ] 


'I  D 0”'D  n-lf 

,01  n.  x-x  x-x 

[x  , x ,...,x  ] 


This  formula  lilies  that  Kergin 's  map  for  the  points  x°,...,xn  is  obtained  by  adding 
to  the  map  corresponding  to  the  points  x°,...,xn  1 the  expression 

/ D 0-”D  n-lf  * 
lx° *"]  X-X 

Sponsored  by  th«  United  States  Any  undtr  Contract  No.  DAAG29-L7&-C-0024. 


I 


This  property  is  reminiscent  of  Newton  divided  difference  expansion  for  polynomial 
interpolation  in  one  variable.  As  we  will  see  this  similarity  is  not  accidental  as  the 
proof  of  (1)  is  based  on  Newton  form  for  polynomial  interpolation. 


Along  with  formula  (1)  we  also  have  available  a convenient  expression  for  the  error 
in  approximating  by  Kargin' s map,  namely. 


(2) 


Pf(x)  - f(x) 


/ 


...  D 


, 0 n.  x-x 

[x,x  ... . ,X  ] 


This  useful  formula  makes  it  possible  to  obtain  error  estimates  for  approximating  by 
Kergin's  map.  For  instance,  since  the  volume  of  s"  is  1/nl  we  have 


I)  Pf-fH 


where  — + — » 1, 
P q 


L (K) 
K ’ 


iCT(VK»n+1 


a 


a*,p 


ll  Dfll 


|a|«n+l 


1/P 


L (K) 


convex  hull  {x^  | 0 <_  j n } and  dq  (K)  - diameter  of  K in  i*. 


Thus  if  F is  a bounded  linear  functional  on  L (K)  which  annihilates  ti  (R  ) then 

n 


|pf  | H Ffl 


(d  (K> ) 


n+1 


( y iDBfi* 

|a|“n+l  L 


(n+W 


1/P 


L (K)  |a|-r.+l  L (K) 

- a quick  ("constructive")  proof  of  the  Bramble -Hilbert  Lemma. 

In  one  dimension  it  is  knotm  that  polynomial  interpolation  to  a fwction  which  is 
analytic  in  a sufficiently  large  region  containing  the  interpolation  points  converges 
geometrically  fast  (as  the  number  of  points  increase)  to  the  function.  Formula  (2) 
provides  us  with  a means  to  extend  this  result  to  higher  dimensions.  He  will  show  that  if 
2 - (x  - (Xj,. ..  ,xfc)  | |x^  | ^1}  (the  unit  cube  in  Rk)  and  f is  a fixiction  on  c"  (Q) 


which  has  an  analytic  extension  to  the  polydisk  D(p)  ” {x  « ( ,...,*„)  1 1 I ± 2p+l) 
p > 1,  then  for  It  c 8 

. k-1 


(3) 


II  Pf  - f H <3 

L (Q> 


(n+k) 


H fit 


D(p) 


(it  seems  to  matter  here  that  Q is  a cube  rather  than  a parallelepiped). 


(n) 


Clearly,  formula  (1)  shows  that  P is  a continuous  map  on  C (K>.  Rergin  proved 


that  when  k < n and  the  points  x°,...,x"  are  in  general  position,  i.e.  every  subset 


of  k+1  points  of  {x°,...,xn)  form  a simplex  of  dimension  k then  P is  a continuous 


map  on  Ck-1(Rk).  This  fact  is  not  obvious  from  (1).  However,  we  will  show  by  elementary 


means  that  when  x ,...,x  are  in  general  position  and  q is  any  differential  operator 


-2- 


* 

.•.a 


'} 


of  order  < n , the  lineer  function*! 


(x°,... ,xnJ 


qf  , ft  cto>  (W 


has  an  extension  as  a bounded  linear  functional  on  C(k_1)  <Rk).  m the  case  Jc-1  such  an 
extension  is  given  by  the  Hermite-Genocchi  formula  for  the  divided  difference  of  a func- 


[ X _ fees#  X ) 
u n 


" [X.tie«|X  ]f  | 


Xj  € R , 


see  Norlund  [4,  p.  16]  (subscripts  will  be  used  for  scalars,  real  or  cr.pl,,).  our  renarks 
concerning  this  aspect  of  Kargin1  s map  leads  us  directly  to  some  simple  properties  of  multi- 
dimensional B-splines  introduced  by  Carl  de  Boor  in  (2].  Needless  to  say,  there  are  many 
intriguing  questions  still  to  be  settled  concerning  these  B-splines. 

Returning  to  formula  (1),  our  first  task  is  to  provide  its  proof.  As  we  said  earlier 

the  Newton  form  for  the  polynomial  of  degree  n which  interpolates  f at  

O n 

(l1)  pf(x)  - f(XQ)  * (x-XQ)  lXQ,X1)f  + ...  + (x-XQ)  . . . (x-Xn_j ) txQ xnlf 

is  the  basis  for  the  proof  of  (1).  it  is  also  well-kno«ei  that  the  remainder  formula  for 

polynomial  interpolation  is  given  by 

(2’)  pf  (x)  - f(x)  - [x,r  ]f  . 

o n 

To  verify  that  Pf  satisfies  the  property  of  Rergin’s  theorem  we  first  observe  that  for  the 
faction  f(x)-e^X  x)  m g(X*x) , into  (1)  gives  Pf  (x)  - p^g(A-x)  where  is  the  polynomial  of 
degree  ^n  which  interpolates  g at  l*x°, . . . , l*xn.  Now,  let  J-  {i/0, . . . .j/1)  **  a 

subset  of  (x  ,...,x  } and  q a homogeneous  differential  operator,  of  order  1 then 


in  i, 
t*  ° * *1 


i'lqUA) 


!*lt) 

i0  a. 

(X*X  , . . . ,X*X  ] 


■ 1 q(iX)[X*x  ,...,X’X  1 )g 
“ i kq(iX)lX*x  ®,...,X»x  *]p  » 


i„ 

[X  °,...,X  lJ 


Since  {e*X  x|x«Rk)  spans  a danse  subset  of  C<n)  <Rk)  we  have  shown  that 


1 


I 


I*  °, . . . ,x  4] 


/ q(f-Pf)  - 0 

i. 


f £CtoV) 


and  thus  P Is  the  Kargin  map. 

The  proof  of  (2)  follows  frosi  (2M  in  a similar  way. 

]( 

For  the  proof  of  (3)  we  require  Cauchy  integral  formula  in  C , see  Bockner  and  Martin 


[1,  p.  33]. 


,,-s  . 2!  f / il£i35 

(2,,i,k  % cjc  (vci,ai+1--(v«k,v 

■ (x1,...,xJ[),  C“  al  - a^  ...  a^l  , 

.1-1  _ 


“l  "k 

*1  ••• 


l°l  “ ^ 


and  i«  tha  contour  {*|  |*|  - 2p*+l}  where  1 < p*  < p.  Thus  we  have 


X « Q , 


and  since 


II D ...  D f«  „ 
X-X°  x-x  L (Q) 


<_  n i x-x  ii  w T ii  D°*fD  m 

IcO  t |o|-n+l  L (0) 

1 ( §->kp~n_l  l al  Ilf*. 

2 |a|-n+l  D(p) 


(3)  follows  from  the  observation  that 

X «»  i l M*  - to*l)l  I 1 1 (dU)  I (n^y^ 

|a|-n+l  | a|«n+l  |a|«n+l  '*  11 

and  that  the  vol  Sn+1  . 1/  (n+1)  I 

in  one  dimension,  polynosiial  interpolation  at  any  (triangular)  sequence  of  points  is 
bound  to  diverge  for  soae  continuous  function  while  for  interpolation  at  the  reros  of  the 
Chebyshev  polynomial  convergence  holds  for  lip  1 functions.  It  would  be  interesting  to 
extend  these  remits  in  Rk  for  Kergln  interpolation. 

We  now  tum  to  the  natter  of  extending  the  linear  ftxictional  defined  in  (4). 

We  begin  by  observing  that 


l 


has  an  extension  as  a bounded  linear  functional  on  C(K)  if  and  only  if  y is  in  the 


The  necessity  of  this  condition  is  a conse 


vector  space  generated  by  x -x 


quence  of  the  following)  Suppose  A e.  R and  A i x-x  , j«l, 
let  f(x)  - ei5*  x « g(6A*x)  then  |Lf|  - |$||yA|  while  II  ftl 


n.  Given  any  6 e R, 


- 1 . Thus  A . 

L~(K) 

If  on  the  other  hand 


then  for  6-1 


t R } spans  a dense  subset  of  C 


Again,  since  {e 


then  L is  a bounded  linear  ftnctional  on  C(K)  for  all 


x are  in  general  position  then  the  above  observation  extends 


by  induction  to  show  that  for  any  y 


is  a sun  of  q the  linear  ftnctional  given  by  (5)  has  a bounded  extension  on  C 


Cne  aay  easily  verify  that,  in  general,  if  q 


(A)  - 0.  Mather  or  not  this 


condition  is  also  sufficient  for  the  boundsAiees  of 


has  not  been  settled 


Finally,  we  conclude  with  some  remarks  on  the  multidimensional  B-splines  introduced 
in  12].  It  is  not  hard  to  see  that  if  x^x0, . . . ,xn-x°  span  Rk  then  there  exists  a 
fmction  M (x  |x°, . . . ,x")  such  that 


. 0 n, 
lx  .... .X  ] 


/f  (x)M(x|x0,. . . ,xn)dx  . 


Note  that  the  expression  on  the  left  is  an  n-fold  integral  over  the  sinplex  S0  while 

0 1 n 

the  right  hand  expression  is  a k-fold  (k  ^ n)  integral  over  the  convex  hull  of  x ,x  . 

Repeated  application  of  (6)  implies  that  if  x^.x1, . . . .x"  are  in  general  position  then 

/ Daf  - / Daf (x)M  (x | x°, . . . ,xn)dx  , |a|  » n-k+1 

[x°,...,xn]  K 

is  a bound  linear  functional  on  C(K).  Moreover,  it  is  a sum  of  integrals  of  f over 
hyperplanes  formed  by  subsets  of  k points  of  x°,...,xn.  Thus  M is  a polynomial  of 
total  degree  <_  n-k  in  every  region  bounded  by  such  hyperplanes  and  globally  has  n-k-1 
continuous  derivatives  on  K,  facts  already  mentioned  in  (2]. 


A formula  for  M(x),  which  reveals  its  geometric  interpretation,  may  be  obtained  by 

"lifting"  the  points  x°,...,xn  into  Rn  . We  let  x°,...,xn  be  any  vectors  in  Rn  such 

that  | - x^  and  that  the  simplex  a determined  by  these  points  is  of  dimension  n. 

R* 

Then 

VOl{u  £ o|u|  “ x} 

,8)  M(x  x°  x»)  . J-  P.7* 5 

,8)  ' n!  vol  a 

n 

(For  k-n,  M (x | x° , . . . , x" ) - ^ /vol  o).  Thus  M is  (with  a different  normalization) 

de  Boor's  multidimensional  B-spline  [2], 

In  spite  of  this  formula's  attractive  geometric  content  it  seems  more  convenient 
to  derive  properties  of  M through  the  defining  equation  (7).  In  particular,  according  to 


formula  (6) 


D M(-|x°,...,xn)  » 1 a. (N(* | 

y 1-1  3 


,x^  1,x^+1, . . . ,xn)-M(*  jx1, . . . ,xn) ) 


y - l a (x3-x°)  . 

1-1  3 

When  k-1,  this  formula  yields  the  equation 


(10) 


M(x|x  ,...,x  ) - M(x|x.  ,...,x  ) 

M'  (x|x x ) - 2 2^1 1 — 


0 n 

(«e®  (21)  and  with  the  initial  condition 


x -x. 
n 0 


•11  the  (univariate)  B-splines  can  be  generated  from  (10).  A similar  remark  applies  for  (9). 

In  the  univariate  case,  an  efficient  and  numerically  stable  aeans  of  generating  B-splines 
is  available, (2).  The  univariate  B-spline  over  its  interval  of  support  is  a nonnegative 
combination  of  two  lower  order  B-splines,  specifically 


(11)  (n-1)  M(x|x  ,...,x  ) - --  M(x|x  ,...,x ) + - — 2.  M (x|x  , . . . ,x  ) 
o n x -x  1 n x -x  l n-1' 

n u n u 

Using  the  identity 


x -*  

(12)  lx„,...,x_]  - **'Jti"'*,Xn-ll+x""-x’"  *r'X0"'"Xn-l^ 


X -X. 


n 0 nO 

for  divided  differences  it  follows  by  the  line  of  reasoning  used  to  prove  (9)  that 

(13)  M(x|x°,x1,. ..  ,xn)  - t M(x|x°,.«.  ,xA  1,s)  ♦ (l-t)M(xlx1,...  ,xn  1,s) 


Z - t x°  + (l-t)xn  . 


In  the  univariate  case,  it  aay  be  easily  proved  that 

0.4)  Mfx|x  ,...,x  ,x)  - — M(x|xn,...,x  ) 

o n ti  0 n 

and  thus  &4)  with  (13)  proves  (11)  by  choosing  x « t Xg  + (l-t)xn>  m the  eulti variate 
case  it  is  highly  plausible  that  a formula  like  (14)  also  holds  and  this  should  lead  to  a 
multidimensional  version  of  (11). 


The  following  theorem  provides  the  details  for  this  intuition, 
n . n 

Theorem  2.  If  x - £ X.x3  , J 1.  ■ 1 

j-0  3 j-0  3 

then  for  n _>  k+1 

n .0  J-l  1+1 


M(x|x°,...,xn)  - -qj-  i l.ll(x|x0,...,x^"1,x3'*■1,...,x,,) 
J-0  3 


Proof.  Me  divide  the  proof  into  several  ancillary  results.  First  we  require  an  observa- 


tion about  divided  differences  which  generalises  formula  (12). 
n n 


Let 


l 3 


j-0 


i V 


I K - 1 

J-0  3 


than 


7- 


(15)  [x0,...,xn)f  - Xj  [*,x0,...,xj_1,Xj+1,...  ,xn)f  . 

Note  that  in  (15)  the  order  of  the  divided  difference  on  each  side  of  the  equal  sign 
is  the  same. 

Proof.  The  right  hand  side  of  (15)  is  a linear  functional  which  agrees  with  the  divided 
difference  appearing  on  the  left  hand  side  on  polynomials  of  degree  n . The  support  of  the 
linear  functional  on  the  right  hand  side  is  contained  in  {xQ, — ,xn,z}.  However,  since  the 
coefficient  of  f(z)  is 

n 1 

Jo  X)  (z-v ■••<*- Vi^'Vi'-'-^-V 


1 p 

- r — r ) X . (z-x . ) - 0 

(Z-XQ)...(Z-Xn)  3 3 

the  support  is  {x„,...,x  } and  thus  the  identity  is  proved. 

On 

We  have  already  demonstrated  in  the  proof  of  (9)  how  identities  for  divided 
may  be  extended  to  higher  dimensions.  This  same  approach  applied  to  (15)  yields 


lx  , 

or  equivalently 


n. 
,x  1 


l 

j'°  [..x°. 


If 


j-1  j+1 

|X  ,x  , 


n, 
»X  ] 


(16) 


M (x|x°, — ,xn)  * l X.M(x|z,x°,... ,x^  1,xj+1,...,xn) 
j-0  3 


differences 
the  identity 


n 4 

z - l X.x3  . 
j-0  3 

The  next  identity  we  require  is 

(17)  M(x|x,x  ,...,x  ) - j^-M(x|x  ,...,x  ). 

We  base  the  proof  of  this  formula  on  the  following  general  consideration. 

Let  K be  a convex  set  Rm  and  y t Rm.  We  form  the  cone  in  Rm  with  base  K x { 0 > 

and  vertex  y * (1),  that  is, 

C - cohull (K  * (0)  u y * {1})  . 


Then 

(18)  vol-+1C  - volnK  . 

We  prove  this  formula  by  first  establishing  it  for  a simplex.  If  K is  the  simplex 
with  vertices  y°,...,ym  then  C is  a simplex  with  vertices  y°  * {o),...,y  * (0),y  * (1). 

The  formula  for  the  volume  of  the  simplex  K is 

-8- 


« 


1 


1 


vol  K “ — 
m ml 


and  applying  this  formula  to  C gives 


vol  C 
m+1 


which  readily  simplifies  to  (18) . 


y y 

m m 


Now,  if  Kj  and  K2  are  disjoint  convex  sets  then  their  corresponding  cones  with  the 
same  vertex  intersect  only  at  the  vertex.  Hence  (18)  extends  to  any  countable  union  of  dis- 
joint simplices  and  in  particular  for  any  convex  set  K . 

0 n-1  n-1  *0  *n-l 

To  prove  (17)  we  first  lift  the  points  x ,...,x  into  R by  means  of  x , . . . ,x 


and  apply  the  relation  (8)  , 


M(x  |x° xn_1)  * 


Voln-k-l{u  £ °IU|  k = X) 
R 

vol  , a 
n-1 


*0  — 1 

where  a is  the  simplex  with  vertices  x , . . . ,x  . The  vector  x is  now  also  lifted  to 
Rn  * in  any  fashion  to  x e R°  * and  then  we  form  the  cone  with  vertex  x x {1}  = x and 
base  o x (o)  and  denote  it  by  o.  It  is  readily  seen  that  the  cone  with  base 
<■  ■ °'V  x)x{0)  and  vertex  x is  {u  € a|u|  ^ ■ x)  and  thus  two  applications  of  formula 
(18)  gives  us 

n , , VOln-k-l{u  £ °lU|k  * X} 

M (x|x x"  ) « — j s— ' « 

1 nl  vol  o 

n 

vol  . {u  e o|u|  . = x} 

_ (n-K)  n~k  _ <Rk 

nl  vol  5 

n 

(n-k)M  ( x | x,x°, . . . ,xn  *)  . 

If  we  now  put  together  formula  (16)  with  (17)  when  x..  = 0,  j = 2 , . . . ,n  and 
. 0 ,1 

x » z » XQx  ♦ Xj^x  we  get 


-9- 


(19)  (n-.k)M(x|x0,. . . ,xn)  * AQM  (xlx1, . . . ,xn)  + X M(x|x°,x2 . ,xn) , 


, 0 ^ , 1 

x ■ XQx  + \^x  . 


With  these  facts  at  hand  we  nay  now  complete  the  proof  of  Theorem  2. 
n 

Suppose  we  write  x - £ X.x' 

j-0  3 

in  th«.  form 


,3 


x = (1  - X )z  + X x , 
n n 

z = X vjxJ  * (i'x«>wj  ■ xj  1 ■ 0,1 n_i  • 

Then  for  every  integer  1 , with  0 <_  1 n-1  we  have 

(n-k)M(x|x,x0,...  ,xl  V1....  ,xn) 

i 0 l-l  l+l  n.  , i 0 l-l  £+1  n-1, 

* (1"^  )M  (X I X 0 • # e $ X |X  f • • • /X  ) + A M (X|  Z|X  t # e • # X »X  # e e # 0 X ) • 

n n 


But  according  to  (16) 
n-1 


r / I 0 1—1  £4-1  nx  „ . i 0 n. 

2,  U M(x|z,x  ,x  ,x  ,...,x  ) - M(x|x  ,...,x  ) 

1*0 


and  hence 


n-1 


(n-k)M(x|x°,...,xn)  « l X M (x|x  , . . . ,x  l,x  ,...,xn) 

£-0 

, "r1  I 0 1-1  1+1  n-1, 

+ X 2 M(XlS,X  • • • • • X |X  t • • • # X ) ■ 

b*0  0 n-1 

Another  application  of  (16)  inplies  that  the  second  sum  above  is  M(x|x  , ...,x  ) and 

thus  the  theorem  is  proved. 

Besides  the  obvious  integral  relation 

M(x|x°,. ..  ,xn)  » i-r-  f e iX  x[iX-x°,. . . ,iX*xn)gdX,  q(t)  « et  , 

(2iri)  Rk 


- 1 / [iX*  (x°-x) i(  X>(xn-x)lgdX 

(2ei)X  Rk 

for  M obtainable  from  (7)  by  setting  f(x)  - eiX*x  and  using  the  Fourier  inversion 
formula  there  is  also  a connection  between  M (x |x°, . . . ,xn) , the  Random  transform  and  the 


univariate  B-spline. 


1 ._k. 


Recall  that  the  Radon  transform  is  a map  of  functions  f e L (R  ) defined  for 
8 t (the  unit  sphere  in  R *)  and  t £ R by 


R0f(t)  - / f(x)dui  (x) 

8»X“t 
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t 


•r  - ■ 


The  Radon  transform  also  has  the  characteristic  property 


I”  <R„f)  <t)g(t)dt  - / f (x)g(6«x)dx,  g«L"(R). 


Thus 


(R.M(-|x°,...,xn))  (t)  - M(t|e  • X0,...,  e *xn) 


a striking  relationship  between  univariate  and  multivariate  B-splines. 
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